Comparison of the density-matrix renormalization group method applied to fractional quantum 

Hall systems in different geometries 
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We report a systematic study of the fractional quantum Hall effect (FQHE) using the density-matrix renor- 
malization group (DMRG) method on two different geometries: the sphere and the cylinder We provide con- 
vergence benchmarks based on model Hamiltonians known to possess exact zero-energy ground states, as well 
as an analysis of the number of sweeps and basis elements that need to be kept in order to achieve the desired 
accuracy. The ground state energies of the Coulomb Hamiltonian ati^— 1/3 and i/ = 5/2 filling are extracted 
and compared with the results obtained by previous DMRG implementations in the literature. A remarkably 
rapid convergence in the cylinder geometry is noted and suggests that this boundary condition is particularly 
suited for the application of the DMRG method to the FQHE. 



I. INTRODUCTION 

Strongly correlated systems in low dimensions are among 
the most active areas in the condensed matter physics. These 
systems contain a large number of particles that interact 
strongly with each other and cannot be understood in a single- 
particle picture. A paradigm of strongly correlated systems is 
the fractional quantum Hall effect' (FQHE) that occurs when 
a system of two-dimensional electrons partially fills one of 
the Landau levels in a strong perpendicular magnetic field. 
Since the kinetic energy is frozen in a partially-filled Lan- 
dau level, the electron-electron interaction is the only rel- 
evant term in the Hamiltonian and leads to the emergence 
of non-perturbative ground-states with fractionalized charge^ 
and anyonic, Abelian and non-Abelian-', statistics. 

Due to the non-perturbative nature of the FQHE, numerical 
methods have played a crucial role since the original work of 
Laughlin^. In particular, exact diagonalization (ED) presented 
itself as a versatile and extremely powerful tool that unraveled 
many of the complexities of FQH systems^. The popularity 
and quick success of ED was due to the specific correlations of 
FQH systems that rapidly minimize the finite-size effects with 
increasing the number of particles in the simulation. Highly 
accurate predictions of the system's properties in the thermo- 
dynamic limit could be obtained by considering systems as 
small as 10 particles^. 

As it is well known, the main bottleneck of ED calculations 
lies in the exponential explosion of the size of the many-body 
Hilbert space as the number of particles grows. While for 
the simplest FQH fractions, such as the Laughlin f = 1/3 
state^, admittedly all essential physical properties can be ob- 
tained in the systems attainable by ED, in the majority of other 
cases ED is not sufficient. This is particularly striking in case 
of spin degree of freedom, or SU(4) internal symmetry if we 
consider FQHE in graphene*'. However, similar constraints 
arise even in the spin-polarized case of the non- Abelian Read- 
Rezayi sequence^, where electrons are believed to pair into 
k > 2-body clusters. Therefore, a non-Abelian iVe-particle 
state at level k is likely to have finite-size effects comparable 
to the Laughlin-like state of N^/k particles. Hence, to address 
the properties of the non-Abelian ground-state it is desirable 



to consider systems at least k times as large. It is therefore of 
essential importance to develop new numerical methods that 
can reach larger system sizes than ED. 

One such method is the density matrix renormalization 
group (DMRG), invented by White^ in 1992. DMRG has 
been quite successful over the last decade when it was ap- 
plied to one-dimensional systems such as the Heisenberg spin 
chains and the one-dimensional Hubbard model. In essence, 
it is a variational method to get the ground state and the low- 
lying energy states of the system. The algorithm contains two 
main parts. One is called the infinite size algorithm which 
grows the system to a big size, and the other one is referred 
to as the finite-size algorithm, which makes the ground state 
converge. The only approximation in the DMRG method is 
the truncation of the Hilbert space according to the eigenval- 
ues of the reduced density matrix for the subblock which is 
used to construct the large system. The more states are kept 
in the reduced density matrix, the higher the accuracy one can 
achieve in principle. It is generally believed that the entan- 
glement entropy of a subregion often grows like the bound- 
ary area of the subregion"^. A larger entanglement entropy, or 
larger correlation, means that one needs to keep more states 
to achieve a sufficient accuracy. The success of DMRG in 
the one-dimensional systems was ensured by the low entan- 
glement between the two subregions which only have a point 
surface between two blocks. 

On the other hand, a FQH system is two-dimensional and 
thus the success of the DMRG method in FQHE is by no 
means obvious. However, in a Landau gauge, the one-body or- 
bitals are Gaussian-localized and provide a mapping to an ef- 
fective "one-dimensional" chain with the long-range Coulomb 
interaction. This motivates an attempt to apply the DMRG 
to the FQHE system. The first such attempt was done for 
the periodic boundary conditions (torus geometry) by Shibata 
et flZ. mostly considering compressible, stripe and bubble, 
phases in higher Landau levels. Feiguin et a/." developed 
a DMRG scheme for the ground state and excited states in 
the spherical geometry for larger systems at filling factors 
ly = 1/3 and f = 5/2. Hard-core interactions were also 
studied on thin cylinders in an unpublished work'^, and for 
bosonic systems in Ref. [3- Most recently, Zhao et al. de- 
veloped an independent DMRG implementation that by far 
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exceeds the previous attempts. In this study, the maximal sys- 
tem size was iVe — 24 for ly — 1/3 and = 34 for i/ = 5/2. 
The independent implementations'^^ appear to differ signif- 
icantly from each other in various aspects, in particular in the 
number of basis states that are kept, ranging from a few hun- 
dred in the torus geometry, up to A^kccp = 20000 states in 
Ref.Sl 

In this paper, we report on the systematic study of the 
FQHE system in the spherical and cylinder geometry based 
on our independent implementation of the DMRG method. 
We address the well-studied FQH systems at fillings = 1/3 
and 1/ = 5/2 with the goal of providing a detailed benchmark 
of the DMRG algorithm and comparing it with the previous 
implementations. New physical results obtained with the cur- 
rent DMRG implementation will be presented elsewhere'^. 

The remainder of this paper is organized as follows. In Sec. 
II we analyze the convergence of the Vi Haldane pseudopo- 
tential"*'^ Hamiltonian on the sphere that is analytically known 
to yield the Laughlin wavefunction as an exact zero-energy 
ground state. In the case of Coulomb interaction, we evalu- 
ate the ground state energy for = 1/3 and ly = 5/2 fillings 
corresponding to the Laughlin^ and Moore-Read^ states. The 
ground-state energies per particle are extrapolated to the ther- 
modynamic limit using finite-size scaling techniques. In Sec. 
Ill, we draw some comparisons with the cylinder geometry, 
which is an alternative geometry for studying the FQHE that 
so far has scarcely been used'*. The convergence for the Vi 
Hamiltonian is found to be significantly faster on the cylinder 
than on the sphere, suggesting that this boundary condition 
might be promising for further studies of the FQHE. Discus- 
sion and conclusions are given in Sec. IV. 



II. SPHERE GEOMETRY 

We study a model for spin-polarized electrons moving on 
the surface of a sphere, with a magnetic monopole 25 placed 
in the center to generate a radially-symmetric magnetic field 
perpendicular to the surface"*. In strong magnetic fields, elec- 
trons in general completely fill (n — 1) single particle Landau 
levels which are considered to be "inert", and all dynamics 
comes from a partially-filled nth Landau level. Any two-body 
Hamiltonian, projected to this nth Landau level (neglecting 
the excitations to higher Landau levels), can be written in the 
usual second-quantized form, 

X! ("Iim2|y|m3m4)a+^a+2am3a,„4. 

mi .m2 .ms ,m4 

(1) 

In the spherical geometry, quantum numbers m/s label the 
z-component of the angular momentum for particle i which 
takes values: — S", ~{S— 1), . . . ,S. The one-body orbitals are 
the monopole harmonics'* Ysim which generalize the usual 
spherical harmonics obtained for 5 = 0. When we target 
a specific many-body state, we also need to adjust the flux 
25 to take into account the so-called shift that determines 
the total number of the available orbitals. This means that 
25= -Ne + S, where 5 is a universal number that character- 
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FIG. 1. (a)The convergence of the ground-state energiesfor the hard- 
core Vi Hamiltonian at = 1/3 as a function of the finite size 
sweeping number when keeping 4000 states in the subsystem. (b)The 
ground-state energy as a function of the number of kept states. We 
perform 10 finite-size sweeps for each point. The energies are on a 
logarithmic scale. 

izes each many-body state, e.g. S = —3 for the Lauglin and 
Moore-Read state. 

Because of rotational and translational invariance, any two- 
body interaction matrix element {m,im2\V\m3m,4) can be de- 
composed asi^ 

{mim2\V\m3m4) = 
2S J 

{Smi,Sm2\JM){Sm3,Sm4\JM)vj^'^/R, 

J=0 AI=-J 

(2) 

where Vj are the Haldane pseudopotentials^ and R = VS£b 
is the radius of the sphere in terms of the magnetic length 
£b — \/ h/eB. The first two terms in the above equation are 
the Clebsch-Gordan coefficients on the sphere. When sym- 
metry is taken into account, at filling = 1/3 the Lanzcos 
method can diagonalize the sparse Hamiltonian matrix for up 
to 14 electrons, corresponding to the Hilbert space dimension 
of - IQS. 

As shown by Haldan©^, the advantage of the pseudopo- 
tential formulation is that model wavefunctions can be de- 
fined as ground states of the truncated Hamiltonians. For 
example, the Laughlin wavefunction is obtained as an exact 
zero-energy ground state for the hard-core interaction with 
Vi > 0, V„iyi = 0, with an excitation gap controlled by the 
magnitude of Vi . From the computational point of view, Vi 
Hamiltonian is nearly as sparse as the full Coulomb Hamilto- 
nian, but it serves as a universal reference to test the accuracy 
of the DMRG code for large systems because the ground-state 
energy is known to be exactly zero for any system size. 

In Fig. [T] we show the ground-state energy convergence for 
different system sizes as a function of the finite-size sweep 
number with the fixed number of states kept (iV^cop = 4000), 
or as a function of the number of kept states in the subsystem 
with a fixed sweep number For a fixed number of kept states. 
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FIG. 2. Ground-state energy per particle for the v = 1/?) state up 
to 20 electrons as a function of size of the system. Solid points 
are the results of exact diagonalization and the blank circles rep- 
resent DMRG results. The two fitting curves are with(upper) and 
without(lower) rescaling the magnetic length as described in Ref. 19 
and|20|.The energy in the thermodynamic limit is ~ — 0.4101e^/^s 
which is consistent with previous studie a"''"* . 



FIG. 4. Finite size scaling of the ground-state energy per particle for 
= 5/2 FQH state up to TVe = 30 electrons. The results are ob- 
tained by keeping up to 5000 states in the truncation and after com- 
pleting 10 finite size sweeps. The energy in the thermodynamic limit 
is consistent before and after rescaling the magnetic length I b , which 
means that the large system sizes have automatically eliminated the 
curvature effects. 
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FIG. 3. The entanglement spectrum for 18 electrons at z/ = 1/3 with 
the Coulomb interaction. The subsystem contains 9 electrons in 25 
orbitals. 



the accuracy of the ground state energy decreases when we in- 
crease the system size. It means that more states for the larger 
systems need to be kept if the same accuracy is demanded. As 
shown in Fig. [lib), increasing the number of kept states ob- 
viously helps the convergence although the energy drops very 
slowly when iVkccp is large. However, for the largest system 
size with = 18 we tested in Fig.[Tl the ground-state energy 
drops to 10^^ when just keeping 4000 states and after finish- 
ing 10 finite size sweeps. This energy is far below the gap 
between the ground state and the first excited state. We as- 
sume the ground state is close enough to the Laughlin state in 
this case. To improve the accuracy, one needs more finite-size 
sweeps and keeping more states in the truncation. 

Having established the convergence scaling for the hard- 
core Hamiltonian, we move to the full Coulomb interaction. 
At I' = 1/3 we calculate the ground-state energy for systems 



up to 20 electrons. All the results are obtained by keeping 
up to 5000 states in the subsystem. With keeping the same 
number of states, we find the efficiency of our code is the 
same as that shown in Fig. 1 of Ref^^. The results for dif- 
ferent system sizes at = 1/3 are summarized in Fig. |2] 
which includes the data both from the ED and the DMRG. 
It shows they match with each other very well. We do the 
finite-size scaling for the ground-state energy per electon with 
a quadratic polynomial, and extrapolate the thermodynamic 
limit energy to be — 0.410048e'^/^B. On finite spheres, it 
has been suggested'''-^" that the curvature effects can be sub- 
stantially minized by rescaling the magnetic length £b- We 
also plot the rescaled energies in Fig. |2] and do the finite-size 
scaling with a linear function. The energy in the thermody- 
namic limit — 0.410136e^/^^ is almost the same as that with- 
out rescaling the magnetic length. This means the large-scale 
study by the DMRG method has already removed the finite- 
size effects coming from the curvature. Our results are also 
consistent with the previous DMRG study" '"*. Besides the 
ground-state energy, we also plot the entanglement spectrum^- 
in Fig. [3] for 18 electrons at ly — 1/3, for which we cut the 
system into two equal parts. The splitting between the confor- 
mal part^' and the non-conformal part, and the counting of the 
conformal states in the entanglement spectrum, demonstrate 
that DMRG has captured the correct topological properties of 
the ground state. 

As a second case, we consider the filling ly = 5/2, believed 
to be described by the Moore-Read Pfaffian stata^. This state 
is more fragile that the Laughlin state and has a smaller gap by 
nearly an order of magnitude. To study the convergence, it is 
in principle possible to use the exact interaction that produces 
the Moore-Read state as a zero-energy ground state, but this 
is much more costly because it is a three-body interaction. 
The results for the Coulomb interaction projected to n = 1 
Landau level are shown in FigH) The ground-state energies 
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FIG. 5. The convergence of the ground-state energy for electrons 
with hard-core interaction on the cylinder, (a) We keep 2000 states 
in the truncation procedure for all the system sizes. The ground-state 
energy drops very fast in the finite-size sweeping procedure, (b) The 
dependence of the ground-state energy on the number of kept states 
for 20 electrons. 



FIG. 6. The mean orbital occupation number for 20 electrons with 
hard-core interaction al i/ — 1/3 on the cylinder. DMRG calcula- 
tion is performed by keeping 3000 states and finishing 6 finite size 
sweeps. 



are obtained for up to 30 electrons by keeping at most 5000 
states. The result for the largest system size presented in this 
plot was obtained within one week on a computer cluster with 
12 cores and 144G memory. With the same scaling techniques 
as in the ly — 1/3 case, we extract the ground-state energy 
per electron in the thermodynamic limit ss —0.3622e^ /£b, 
consistent with Ref. 1 liandJ4i. 



III. CYLINDER GEOMETRY 

To complement the results obtained in the spherical ge- 
ometry, in this Section we consider the cylinder geometryi^. 
Cylinder geometry is interesting because it shares some fea- 
tures with the compact geometries, such as sphere or torus, 
but also possesses two open boundaries, which makes it con- 
venient for the study of the edge effects, like the disk geome- 
j.j.jj22-24 Compared to the sphere, the attractive feature of the 
cylinder is the flat surface and lack of curvature effects. 

Cylinder boundary condition is compatible with the Landau 
gauge where periodic boundary condition in assumed along 
one direction (say y-axis) with a repeat distance L, and open 
boundary condition in the other direction (x-axis). The single- 
body wavefunction in the lowest Landau level is given by 

i;^{x,y) = , ^^^y^-i.^^..^sfli(^^^ (3) 

where fc„i = l-KmjL is the momentum for the mth orbital. 
The orbital index to takes values 0,1,..., Norh ~ 1, and the 
distance in x-direction between two nearest orbitals is 2tt/ L. 

For a finite size system with particles at filling v — l/Z 
for example, the number of orbitals is Norb — 3A^e — 2, and 
thus the area of the cylinder is quantized to be 2TTNorb^%- To 
accommodate the finite number of the orbitals Norb, we fix 
the extent in the x-direction to be X ~ 2TiNorb£% I ^- Similar 



to the torus geometry, the properties of many-body states de- 
pend on the aspect ratio A = X/L — 2TTNorbi%/L'^- In the 
following we concentrate on the FQH states that are realized 
in the vicinity of A = 1. 

For the Vi hard-core interaction, the Hamiltonian (in units 
£b ~ 1) can be written in a simple formi^: 

^ = ^ E - n2)e-("^+"^)/2a+_,,a+^,a,„+„+,a^.4) 

m,n,/ 

Exact diagonalization studie a'^''^ show that H has a zero- 
energy ground state, but the nature of the ground state changes 
from the incompressible liquid to the charge-density wave 
upon varying the aspect ratio. Here we focus on the liquid 
state and obtain it by DMRG method for large systems up to 
20 electrons with the hard-core interaction. 

The convergence of the ground state energy as a function 
of the sweep number and the number of kept states is shown 
in Fig. |5] For the system with 20 electrons, we obtain the 
ground state energy 10^^'^ when keeping only 2000 states and 
after completing 6 finite-size sweeps. On the other hand, if 
we look at the final ground state energy as a function of the 
number of the kept states, we observe that the same accuracy 
can be reached even by keeping only 1000 states. 

To verify that the ground state is indeed the Laughlin state, 
we plot the average occupation number (c+ c.,„) for the system 
with 20 electrons in 58 orbitals in Fig.|6l For an incompress- 
ible liquid, the average occupation number is roughly constant 
in the bulk and equal to with some deviations close to the 
two edges. This is indeed what we observe in Fig.|6] 

We also plot the entanglement spectrum on the cylinder. 
Fig- 111 Because we use the ground-state of the Vi pseudopo- 
tential Hamiltonian, the entanglement spectrum only contains 
a conformal branch, but is otherwise similar to the spectrum 
obtained on the sphere. Fig. [3] In particular, the counting of 
the conformal levels is identical in the two cases (up to the 
limit set by the size of the sphere). Note that although the 
true energy spectrum reflects the presence of two edges on the 
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FIG. 7. The entanglement spectrum of the Laughlin state for 20 elec- 
trons on the cylinder Because of the hard-core interaction, the en- 
tanglement spectrum only contains the conformal branch, with the 
same counting as in Fig. [3] 



cylinder, the entanglement spectrum involves only a single cut 
and thus probes only a single edge, in complete analogy with 
the sphere. 



tion of the DMRG method. In the spherical geometry, the 
DMRG results for the ground state energies at filling v —Xji 
and V = 5/2 are consistent with the exact diagonalization 
study for small system sizes, and the previous DMRG stud- 
jg^iiji4 j^-gg system sizes. For the largest system size we 
have reached, the eiTor of the ground state energy is about 
10^'* which is roughly two-three orders of magnitude below 
the energy gap to the excited states. The consistency in the ex- 
trapolation of the ground-state energy shows that these system 
sizes have negligible curvature effects. 

The application of the DMRG method to the cylinder ge- 
ometry shows much higher efficiency compared to the sphere. 
Based on the convergence for the Vi interaction, we expect the 
cylinder to be the more promising venue for the future appli- 
cations of DMRG. Due to the presence of two open edges, the 
treatment of the full Coulomb interaction is not as straight- 
forward as in the compact spherical geometry, and requires 
special care in defining the confining potential to contain the 
fluid. One may furthermore expect various phase transitions 
as a function of the aspect ratio and the magnitude of the con- 
fining potential relative to e^eis- Details of these studies will 
be presented elsewherai^. 



IV. CONCLUSIONS AND DISCUSSION 

We presented a systematic study of the FQHE at two well- 
known and important filling factors, = 1/3 and v = 5/2, 
for boundary conditions using our independent implementa- 
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